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j^ ■ Abstract 

We use a high-precision Monte Carlo simulation to determine the universal 
specific- heat amplitude ratio A^/A_ in the three-dimensional Ising model via the 
impact angle (j) of complex temperature zeros. We also measure the correlation- 
length critical exponent u from finite-size scaling, and the specific-heat exponent 
a through hyperscaling. Extrapolations to the thermodynamic limit yield (j) = 
59.2(1.0)°, A+/A_ = 0.56(3), u = 0.63048(32) and a = 0.1086(10). These results 
are compatible with some previous estimates from a variety of sources and rule out 
recently conjectured exact values. 



1 Introduction 

Although it is one of the most investigated models of statistical mechanics, the still- 
unsolved Ising model in three dimensions has recently again come under the spotlight 
[IlElEllllElinilTlElElIinillllllS]. This is because of (i) controversial claims as to 
the exact values of the critical exponents from analytical means and (ii) advances in 
algorithmic approaches which have allowed for greater numerical precision. 

On the analytic side, two relatively recent papers claim to have found the exact expo- 
nents of the 3D Ising model [H [2]. They claim rational values for the critical exponents, 
including a = and u = 2/3 [5]. A speculative formula embracing the results for the 
magnetic exponent 6 was subsequently given in Ref.[l] and the conjectured values for crit- 
ical parameters in three dimensions were used in Ref.[5]. Despite criticism in Refs.[6l [7], 
experimental support for these critical exponents has been claimed in Ref. [8]. Further 
theoretical arguments in favour of these exponents are given in Ref. [5] (see also Ref.|lU]). 
Refs. [H [21 E] also advocate the existence of a multiplicative logarithmic correction in the 
critical behaviour of the specific heat in the model. 

Given the intriguing coincidence of the claimed values of the critical exponents in 
Ref.[T] and Ref. [2], and their connection with the Rosengren conjecture [13], precise mea- 
surement of a and u are of renewed interest. Precision is important here, not least because 
if any disagreement with RG is found, this will have serious implications for one of the 
foundation stones of theoretical physics. This is also why the claims of Refs. [H [2] caused 
so much controversy - they are not in agreement with precise, though non-exact, results 
from other methods. Moreover, while the values of the critical exponents claimed in 
Refs.[T], [2] obey the standard scaling relations, a logarithmic term in the specific heat 
would violate the scaling relations for logarithmic corrections unless certain conditions 
hold for the for the locus of Fisher zeros [13] . 

Although critical exponents have more commonly been used, amplitude ratios of- 
ten offer a more discerning way to identify universality classes [15]. These have again 
recently received increased attention in the case of the of the 3D Ising model because 
improved numerical techniques now permit more stringent investigations. The most re- 
cent Monte Carlo simulations have yielded the estimates A^/A^ = 0.532(7) [TT] and 
Ajf^/A^ = 0.536(2) [12] and these compare well with previous estimates from a variety 
of approaches in Table 1. The first entry in the table is conspicuous in that (i) it is 
significantly different to the other estimates and (ii) it comes from measuring the angle 
of impact of complex-temperature Fisher zeros onto the critical point, instead of through 
the more direct approaches employed by the other references cited. Therefore we con- 
sidered it worthwhile revisiting this method, making use of the greater computational 
power available today, in an effort to determine whether the inconsistency is indicative 
of a deeper problem or whether it is simply down to limited numerics. For these reasons, 
our analysis focuses on the scaling and impact angles. 

Our detailed numerical analysis upholds previous estimates for the critical exponents 
coming from a variety of approaches. It also brings specific-heat amplitude-ratio estimates 
from the impact-angle approach into line with other estimates and rules out the recent 
conjectures for exact exponents. 



2 Model and observables 

The partition function of the pure 3D Ising model on a lattice of extent L without an 
applied magnetic field and at an inverse temperature /3 = 1/T is 



^i(/3) = X] exp /3 ^ GiGj 
W,} \ {id) 



(2.1) 



where the sum over configurations {ai} is taken over Ising spins, ai G {±1}, and where 
{i,j) denotes nearest neighbours. This can also be expressed in terms of the total energy, 

E, as 



ZL{/3) = J2p{E,P)t 



-PE 



f2.2) 



where p{E, (3) is the density of states. We define the reduced temperature as a dimen- 
sionless distance from criticality, 

«=1-|. (2.3) 

where /3c is the critical values of /3. The specific heat of the system, is defined as 



C=ym)p-{E)l), 



(2.4) 



where V = Lr" is the volume of a cubic system of linear extent L, and (. . .)^ denotes the 
thermal average at /3. Close to a continuous phase transition, the specific heat scales with 
reduced temperature t as 

C{t)c^A^\t\-\ (2.5) 

Table 1: Measurements of the amplitude ratio for the specific heat of the Ising model in 
three dimensions. The asterisk indicates that Marinari's Monte Carlo estimate was from 
the impact angle of Fisher zeros. 



Authors & Reference 


Year 


Method 


^+M- 


Marinari [16] 


1984 


Monte Carlo* 


0.45(7) 


Belanger and Yoshizawa |17 


1987 


Experiment 


0.53(1) 


Bagnuls et al |18j 


1987 


Field theory 


0.541(14) 


Liu & Fisher [I9] 


1989 


Series expansions 


0.523(9) 


Hasenbusch & Finn |20j 


1998 


Monte Carlo 


0.560(10) 


JNowicki et al |21j 


2001 


Experiment 


0.536(5) 


Campostrini et al |22j 


2002 


Series expansions 


0.532(3) 


Canfora et al 23] 


2009 


Phenomonology 


0.54 


Feng & Blote [H 


2010 


Monte Carlo 


0.532(7) 


Hasenbusch 12 


2010 


Monte Carlo 


0.536(2) 



where the amphtudes /!+ and y4_ correspond to t > and t < respectively. The 
susceptibihty is also of interest and can be defined as 



yKV^J- /P - V^-^lp) 



X = -{{M^)p-{M)l), (2.6) 



where M = Y2^ ai is the magnetisation. One can estimate, for instance, the critical point 
using the scaling of the specific heat maxima with the lattice size. The currently widely 
accepted value [24] for this point is f3c{L = oo) = 0.2216546. 

A complex zero in the partition function indicates a non-analyticity in the free energy. 
In the thermodynamic limit {L -^ oo), the pinching of such zeros of the real temperature 
axis precipitates a phase transition. Nevertheless we can study the transition in a finite 
system by analytic continuation to complex temperatures, /3 = r] + i^. In this case the 
partition function includes both oscillating and damping factors: 

Z(/3) = 5^p(E,/3)e-(''+^«)^ = J2piE,/3)e-'^''[cos{^E)-zsm{^E)] . (2.7) 

E E 

Rescaling with Z[Re(/3)] we define the quantity: 

T^(r.^\ - ^(/^) - Y.EPiE,P)e-^^[cos{iE)-is\n{iE)] 
""^"^'^^ - Z[Re(/3)]- Y.EPiE^P)e-^^ ^''"^ 

= (cos(e^)),-^(sin(e^)),. (2.9) 

Therefore the partition function at complex temperature P = rj + iC, can be constructed 
using expectation values taken at real temperatures P = rj. A Fisher zero is a complex 
temperature value [r], ^) such that R = 0. These values can be ordered by their distance 
from (3c{oo) and we denote the first and second zeros by 13^^^ and (3^'^\ respectively. 

Our objective is to obtain the locations of the Fisher zeros for several system sizes. 
To achieve this we performed MC simulations in the vicinity of the critical temperature. 
The crucial point is that we can use histogram reweighting techniques [25] to evaluate 
{cos{^E))rj and (sin(^£'))^ near the simulation temperature. With this in hand, we can 
then minimise \R\ using numerical optimisation methods [26] . 

3 Finite-Size Scaling of the Fisher Zeros 

In Ref.[27] the finite-size scaling (FSS) of the complex thermal coupling (or Fisher) zeros 
was obtained. The zeros nearest to the real axis (3^^^ scale with lattice size as 

/?(!) = L-i/V'i(0), (3.1) 

where /~^(0) is a complex number in general. So, both the real and imaginary parts of 
the Fisher zeros should scale with the same exponent l/u. In Appendix A, we discuss an 
exception of this behavior, and determine under what conditions the real and imaginary 
parts scale with different power laws. 



In addition there is a relationship between the impact angle of Fisher zeros and the 
critical amplitudes of the specific heat. If is the angle between the negative sense of the 
real temperature-parameter axis and the locus of nearby Fisher zeros in the thermody- 
namic limit, this relationship is 

, ,, cosfvra;) — y4_M_i_ 
tan 2 - a = ^ ' / + . 3.2 

Notice that if, instead, we define the impact angle as between the positive sense of the 
real temperature axis and the locus of zeros then we obtain a similar equation with the 
roles of A_ and A^ interchanged: 

tan 2-aU= ^ — '- — ^ . 3.3 

sin(7ra) 

It is clear from either equation that, since the critical exponent a and the amplitude 
ratio A_/y4+ are universal, the impact angle is also an universal quantity. If = 7r/2, 
these equations give that A+ = A_, which is physically reasonable; vertical impact signals 
symmetry on either side of the phase transition and equality of amplitudes. 

In this paper, we will follow the first parametrization of the impact angle: > will be 
the angle between the negative direction of the real temperature axis and the Fisher-zeros 
locus. We will then use the relationship (13. 2 p to determine the amplitude ratio A_/A^ in 
a manner independent to [16], and different from Refs. [HI [121 [TTl [HI [IS |20l [2ll |22] (see 
Table 1). We will show that the different methods yield compatible results. 

4 Simulation Details 

We performed extensive simulations of the model for linear lattice sizes from L = 4 to 
L = 72 with periodic boundary conditions. The spins were updated by combining ten 
Wolff single-cluster algorithms with a Metropolis one. We define in this way our unitary 
Monte Carlo sweep (MCS). We measured after every MCS, performing 10'^ measures 
for each lattice size after thermalisation. We used a 64-bits pseudo-random generator 
combining a Parisi-Rapuano wheel with a congruential one [28] . 

We simulated twenty independent pseudo-samples for each lattice size starting from 
random configurations. We merged their MC histories checking that every pseudo-sample 
was thermalised. To check this point we can plot the logarithmic binning of some quanti- 
ties for different samples and see the clear compatibility of the plateaus between samples. 
This is illustrated in Fig. [Tjfor our largest system, for both the specific heat (12. 4p and the 
magnetic susceptibility (12. 6p . 

We have performed two set of simulations. Firstly we simulated every lattice size at 
(3c{L = oo) = 0.2216546 [j, as this should be the more direct approach to estimate critical 
quantities. Nevertheless, we observed unusually large error bars (always computed via 
jackknife blocking between samples) for the locations of the second zeros in the larger 



'^In terms of the associated variable u = exp(— 4/3), this is given by Uc = exp(— 4/3c) = 0.4120468. 
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Figure 1: (Color online) Log-binning of susceptibility and specific heat for five random 
pseudosamples with L = 72. Error bars are typical deviations inside each bin. The first 
block only includes seven measurements, explaining the deviations for small times. 



lattices. We then performed simulations closer to the previously estimated locations of 
these second zeros for every L > 24. We call these new temperature set {f3!,[^}. We 
obtained fully compatible values for the first zero location from these new temperatures 



Table 2: Simulation details for each system size L. Nth denotes the number of MCS (in 
units of 10^) performed in each pseudosample before taking measurements. We simulated 
every L at /3c(oo) = 0.2216546. In addition, in the third column, we present the second 
simulation temperature used for the largest lattices, (3i.^^. 



L 


iVth 




4 




- 


6 




- 


8 




- 


12 




- 


16 




- 


20 


5 


- 


24 


5 


0.22390 


32 


10 


0.22300 


48 


10 


0.22259 


56 


30 


0.22239 


64 


40 


0.22213 


72 


50 


0.22200 



0.01 




0.215 0.217 0.219 0.221 0.223 0.225 0.227 
ry = Re(/3) 

Figure 2: (Color online) Graphical analysis for L = 16 with 2 x 10^ measurements and 
a 1000 X 1000-point extrapolation grid. Dark (red) lines represent the points where 
{cos{C,E))tj changes sign, light (green) lines are the corresponding signals for (sin(^£'))^. 
A crossing of the lines is an imaginary zero of the partition function. The extrapolation 
range allowed the location of the first two zeros. 



and we strongly reduced the error bars for the second zeros. It provided us a clear evidence 
of the sensibility of the second zero location and discouraged us from studying higher- 
index zeros. In Table |2] we give the concrete simulation information for each system size. 
The total estimated simulation time is around a year of a 3GHz CPU. 



5 Results 

To determine the values of the universal amplitude ratios we first obtained the Fisher 
zeros which lie closest to the real temperature axis for each system size. We used the 
two-step procedure described in Ref.[2H], that is: graphical estimation of the approximate 
zero location followed by a multidimensional numerical minimization starting from the 
previously estimated point. 

For the graphical estimation we had to perform an extrapolation in Re(/3) = t] (via 
histogram reweighting, see for example, Ref. [22]) and Im(/3) = ^. Our goal is to obtain 
the points where (cos(^-E'))^ and (sin(^£'))^ vanish simultaneously. We defined the secure 
T) extrapolation range using the square-root of the energy fluctuations as showed in the 
reweighting method [29]. For the optimum C, range we used the known scaling of the 
imaginary parts of the zeros with the system size. A compromise must be achieved between 
resolution (number of extrapolated points) and computation time. The analysis time for 
a lattice with 2 x 10^ measures was around a week of a 3GHz CPU for a extrapolation 
grid with 300 x 300 points. In Fig. [2] we show the output of one of our analyses. 



Table 3: First and second partition function zeros in terms oi u = exp (—4/3), where /3 is 
the inverse temperature. Notice that if u is a zero its complex conjugate, u*, is also. 



L 


Re{u^^^) 


lm{u^^^) 


Re{u^''^) 


lm{u^''^) 


4 


0.3842660(79) 


-0.0877447(77) 


0.344594(055) 


-0.143306(049) 


6 


0.3975600(88) 


-0.0454084(75) 


0.377482(124) 


-0.072941(165) 


8 


0.4027230(57) 


-0.0285868(44) 


0.390136(159) 


-0.045566(171) 


12 


0.4070150(32) 


-0.0149382(25) 


0.400370(224) 


-0.023291(188) 


16 


0.4088090(27) 


-0.0094326(32) 


0.404832(156) 


-0.014731(191) 


20 


0.4097530(30) 


-0.0066037(20) 


0.406870(211) 


-0.010347(227) 


24 


0.4103170(18) 


-0.0049403(18) 


0.408153(030) 


-0.007684(030) 


32 


0.4109431(15) 


-0.0031260(20) 


0.409582(034) 


-0.004871(034) 


48 


0.4114612(09) 


-0.0016399(08) 


0.410759(037) 


-0.002597(038) 


56 


0.4115855(06) 


-0.0012821(10) 


0.411005(022) 


-0.001995(026) 


64 


0.4116732(07) 


-0.0010367(07) 


0.411201(015) 


-0.001615(022) 


72 


0.4117365(07) 


-0.0008596(06) 


0.411343(016) 


-0.001324(014) 



Once we knew the approximate locus of a given zero, we performed an AMOEBA 
multidimensional minimization [30] starting from this point. In this way we computed 
the zeros location with quite high accuracy. In Table [3] we display the results for the 
locations of the first two zeros for each lattice size in terms of the conjugate variable 
u = exp (—4/3). All error bars are computed defining jackknife blocks between the twenty 
pseudo-samples. The larger error bars for the estimations of the second zeros for the 
smallest lattices (L < 20) are due to the longer temperature extrapolation, see table O 
In order to define the impact angle we use the locus with Im/3 > or Irau < 0. 



5.1 Critical Exponents 

As was stated in Secj3l we can obtain the critical exponent u from the scaling of the 
real and imaginary part of the zeros (see Appendix). The scaling of the first two zeros 
is shown on a double-logarithmic scale in Fig. 121 As an additional estimator, we can use 
the distance from the complex zero to the critical point, \u — ud- Therefore we fit each 
of Im(M(-'^), Uc — Re(M*^-^^) and \u^^^ — Uc\ to the form 



aL 



-\lv. 



1 + bL- 



(5.1) 



which includes the leading correction-to-scaling term. Firstly we tried to obtain both u 
and u from fits of our data sets. We were not able to obtain them simultaneously with 
the same high accuracy as that in Ref. [31]. Eliminating u by taking ratios and fitting 



Im(u(^))/|Re(u(^)) - ud to A + BL''^ gives the estimate oo = 0.77(9) using L = 12 to 
L = 72 or u = 0.63(16) using L = 16 to L = 72. These compare reasonably to the 
estimate u = 0.832(6) in Ref. [3T] . Next we decided to focus on the computation of the 
critical exponent u. To this end, we fixed the value u = 0.832(6) [31] and obtained the 



Table 4: Estimated values for the critical exponent v from different estimators (first 
column). The first error bar correspond to the statistical error, while the second one is 
due to the uncertainty in uj. We give in brackets the size range used in each fit. 





«(1) 


u^') 


Re(M) — Uc 
Im(M) 

U — U(\ 


[8-72]: 0.6328(12) (2) 
[6-72]: 0.63048(25) (7) 
[8-72]: 0.63052(41) (2) 


[4-72]: 0.6576(82) (5) 
[4-72]: 0.6344(44) (3) 
[4-72]: 0.6388(39) (3) 



estimates for v shown in Table HI We discarded data from small sizes until we reached 
a good value for the confidence level (CL) of the fito To further ensure our control of 
the corrections to scaling, we checked that our fits are compatible with the fits obtained 
discarding the next smallest L value. We give the error bar corresponding to these latter 
fits. 

Taking our best estimation (the one with the smaller error bar), we obtain the value |j: 



0.63048(32) 



(5.2) 



This is in very good agreement with the world mean value provided in Ref. [23] (namely 
V = 0.6301(4)). For comparison, the estimate coming from Ref. [31] is i^ = 0.63002(10). 

The specific-heat critical exponent can then be estimated using the hyperscaling rela- 
tion 

ud = 2-a, (5.3) 

where d = 3 in the present case. We find 



a 



0.1086(10). 



IS a 



(5.4) 
0.1097(12) and the estimate from 



The world average coming from Ref. 
from Ref. [31] is 0.10994(30). 

The conjectures regarding exact values in Refs.[Tl|2] include a logarithmic divergence 
in the specific heat. According to the scaling relations for logarithmic corrections, such a 
logarithm would necessitate the vanishing of a and an impact angle other than 45° [14] . 
Since the estimate for a in Eq. (l5.4p is over 100 standard deviations away from zero, these 
conjecture are unsupported. 



5.2 Impact Angle and Amplitude Ratios 

We also considered it interesting to contrast the location of the Fisher zeros with the 
temperatures where the specific heat is maximum. For every lattice size, we extrapolated 



^The CL is the probability that x^ would be bigger than the observed value, supposing that the 
statistical model is correct. As a rule, we consider a fit insufficient whenever CL< 10%. 

^ Another test of the robustness of our data is to perform the weighted mean of the data in Table H] 
(where the weights are their inverse variances). In this way we obtain 77 = 0.63069(29) and a = 0.1079(9). 



8 




Figure 3: (Color online) Top: Scaling of the distance from Uc(oo) = 0.41204677 of the 
real part of the first (blue circles) and second (red squares) Fisher zeros. Bottom: Scaling 
of the imaginary part of the first (blue circles) and second (red squares) Fisher zeros. In 
each case the error bars are smaller than the symbols. 



Eq. 02. 4p from the simulated temperature to obtain its maximum, see Fig. |H In Fig. |5]we 
plot the locations of these maxima and the locations of the Fisher zeros for every lattice 
size. From this plot one can see that the angle between the positive sense of the real axis, 
the specific-heat pseudocritical point and the first Fisher zero change significantly as the 
lattice size increases. 

There are very few works in the literature which use the impact angle of the Fisher 
zeros to determine amplitude ratios. Those which do exist are rather old and only had 
access to very small lattice sizes. Marinari's estimate [16] A^/A_ = 0.45(7) is based on 
lattices of extent L = 4 to L = 8. An additional motivation for our investigations is to 
re-examine the impact angle in the light of the much larger lattices available to us today. 

Unfortunately, there is no FSS theory for the impact angle for Fisher zeroqj. However, 



One could attempt to build a FSS theory for the impact angle by noting that, in the infinite- volume 




0.21 0.215 0.22 0.225 0.23 0.235 0.24 0.245 0.25 

H 

Figure 4: (Color online) Specific heat as a function of the inverse temperature for different 
system sizes. We can see the clear approach to the asymptotic value /3c(oo) = 0.2216546 
as the system size increases. 



Table 5: Angles formed between the first and second zeros and the real axis in the u and 
/3-planes. 



L 


01,c(m) 


01,c(/3) 


02,c(m) 


02,c(/3) 


01,2(m) 


01,2(/3) 


4 


72.431(5) 


78.815(5) 


64.79(2) 


75.90(2) 


54.47(5) 


72.15(5) 


6 


72.31(1) 


75.55(1) 


64.64(9) 


70.1(1) 


53.9(2) 


62.6(2) 


8 


71.93(1) 


73.96(1) 


64.3(2) 


67.6(2) 


53.4(5) 


58.8(5) 


12 


71.38(1) 


72.43(1) 


63.3(5) 


65.0(5) 


51.5(9) 


54.2(9) 


16 


71.05(2) 


71.72(2) 


63.9(6) 


65.0(6) 


53(2) 


55(2) 


20 


70.84(2) 


71.30(2) 


63(1) 


64(1) 


52(2) 


54(2) 


24 


70.72(2) 


71.06(2) 


63.1(2) 


63.7(2) 


51.7(4) 


53.6(4) 


32 


70.54(3) 


70.76(3) 


63.2(4) 


63.5(3) 


52.1(8) 


52.6(8) 


48 


70.34(6) 


70.45(6) 


63.6(7) 


63.8(7) 


54(2) 


54(2) 


56 


70.38(8) 


70.47(9) 


62.4(6) 


62.6(6) 


51(2) 


51(2) 


64 


70.16(6) 


70.24(5) 


62.4(5) 


62.5(5) 


51(1) 


51(1) 


72 


70.13(4) 


70.19(3) 


62.0(6) 


62.1(6) 


50(2) 


50(2) 



a number of ways to estimate the impact angle suggest themselves. These include the 



limit, one may describe the locus of zeros as ti(r) = u^+rcxp (i0(r)). Taylor expanding 0(r) = (p+rtj)' - 
one finds that the impact angle is given by 



tani 



Imw(r) 
Rcu(r) — Uc 



rcj)' secret) + 0{r'^). 



(5.5) 
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Figure 5: (Color online) First two zeros as a function of the system size. The locations 
of the pseudocritical temperatures from the maximum of the specific heat are also given 
on real axis. We indicate with a large red disc the asymptotic value Mc(oo). The inset is 
a magnification for larger lattices. 



angle at the critical point between the negative sense of the real axis and the first or 
second Fisher zero as well as the angle formed by the impact of the line joining these 
zeros and the real axis. In Table [5], each of these are listed in both the u- and /3-planes. 
The notation is as follows: 



[u 



<yj,c\"'j — angle at Uc between the j zero and the real w-axis. 

• 01,2 (w) = angle between the line joining u^^\L) and u^'^^L) and the real w-axis. 

• (pj^cW) = angle at /3c between the j*^ zero and the real /3-axis. 

• 01,2 (/3) = angle between the line joining (3^^\L) and (3^'^\L) and the real /3-axis. 

Angles measured in the u- and /3- planes should converge since the transformation from /3 
to M is a conformal one, which preserves angles locally. However, for small lattices, where 
the zeros are relatively far apart, we may expect different impact angles for the /3- and 
M-planes. But in the large L region and when the zeros are sufficiently close to the real 
axis, the conformality of the transformation implies that both angles should be nearly 
equal. This fact serves as a check that we are approaching the asymptotic region in the 



To achieve the known FSS of the zeros u'^^(L) [27], write rj{L) ~ R [j /L 

tan (j) = tan (j)j (L) — pjL^~ + . . . , 



where (/>j(L) = Imu^^^L) / {Reu^^'> (L) - Uc) and u^^^L) = u{rj{L)) and pj = Rj^/'"^4)' sec 
nately we have found this approach fails due to strong corrections. 



Replacing r by rj{L), 

(5.6) 
Unfortu- 
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computation of the impact angle and we can safely discard the results from lattices that 
are too small to achieve approximate conformality. We find that approximate conformality 
kicks in at around L > 32. This means that the earlier measurements of Refs.|16[[27] are 
unreliable and justifies our re- visitation of the problem. The results for the various angles 
and the various lattice sizes are given in Table El The angles were obtained for the data 
simulated at {f^i.^^} to use the correlations between the two zeros to reduce the errors via 
the jackknife method. 

In the absence of a FSS theory for impact angles, we conservatively take the wedge 
formed by the smallest and largest angles for which the lattices exhibit approximate 
conformality to approximate the approach of the zeros to the real axis. This wedge can 
be used to estimate upper and lower bounds for the impact of the locus of Fisher zeros 
onto the real axis. The impact angle may then be approximated by the average of these 
bounds, with their difference delivering an error estimate. In Table [6l these averages are 
listed for lattice sizes that exhibit the desired approximate conformality. 

From the values listed in Table El we observe a drift towards smaller values of as the 
lattice size increases. The values for the largest lattice studied (L = 72) and its associated 
estimate for the amplitude ratio are 

= 59.9(8)°, ^ = 0.58(2). (5.7) 

On may alternatively attempt to extrapolate to the thermodynamic limit using a correc- 
tion ansat^j 

0(L) = + 6L-'^ , (5.8) 

where 4>{L) is the impact angle computed with the data from the lattice size L and is 
the asymptotic value of this angle. 

With u = 0.832(6) from (^, we obtain (fitting only L > 48) 

= 59.2(1.0)°. (5.9) 

The corresponding estimate for the amplitude ratio is 

^ 0.56(3) . (5.10) 



A. 



This compares with Hasenbusch's estimate A^/A^ = 0.536(2) from Ref. [12] (see also 
Table [ID . 



^The rationale for this fit is the following. We know that 

tan0(L) = (fm(u(2)) _ l-ai{u^^^))/{Rc{u''^^) - Re{u^^^)) . 

The scaling of the numerator and denominator is given by eq. (|5.1|) . and finally assuming small angles 
we obtain 

= 0(L) + O(L-"). 
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6 Conclusions 

Because of recent claims around the values of the critical exponents of the three-dimensional 
Ising model, and in the light of renewed interest in the measurements of amplitude ratios 
in this model, we decided to investigate this well-studied model from an alternative per- 
spective, namely that of complex-temperature zeros. In particular, the manner in which 
such zeros impact onto the real axis dictates the amplitude ratios for the specific heat. 
Only very old (almost 3 decades old) measurements of this impact angle exist in the litera- 
ture, so another motivation was to investigate this using the powerful algorithms and large 
lattices currently at our disposal. Indeed, we have found that previous measurements of 
the impact angle were made on lattices which were too small to manifest conformality in 
the complex temperature plane. 

In the absence of a finite-size scaling theory for the impact angle, we were able to 
estimate upper and lower bounds for its value. These lead to an estimate (f) = 59.2(1.0)°, 
which in turn leads to an estimate for the amplitude ratio of y4_|_/y4_0.56(3). In addi- 
tion to the amplitude ratios, we investigated the correlation-length critical exponent and 
determined that u = 0.63048(32). The corresponding value for the specific-heat critical 
exponent is a = 0.1086(10), from hyperscaling. These values are compatible with the 
most accurate recent measurements, and therefore reinforce confidence in our estimates 
for the impact angle. All of our estimates are consistent with the values currently widely 
accepted and are incompatible with the values recently claimed in Refs. [H [2] . 
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acknowlege support from Research Contracts No. FIS2007-60977 (MICINN), GR10158 
(Junta de Extremadura) and ACCVII-08 (UEX). 



Table 6: Impact angle estimation from the average between 0i,c('u) and 0i,2('u) only for 
the sizes respecting conformality. We also estimate the universal amplitude ratio Aj^/A- 
from this angle. 



L 


Angle Average 


\.^+ / ^— ) average 


32 


61.3(4) 


0.610(9) 


48 


62(1) 


0.63(2) 


56 


60.5(8) 


0.59(2) 


64 


60.5(7) 


0.59(2) 


72 


59.9(8) 


0.58(2) 
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A Appendix 

In FSS theory, the shift exponent A characterizes the scahng of the pseudocritical point 
towards the critical point as the thermodynamic hmit is approached. The pseudocritical 
point is given by the location of the specific-heat maximum or by the real part of the 
first Fisher zero and the approach is as L~^. In most models exhibiting higher-order 
phase transitions, the shift exponent A coincides with the inverse of the correlation-length 
critical exponent l/i^, but this is not a direct conclusion of FSS theory and is not always 
true. 

For example, for the Ising model in two dimensions, Ferdinand and Fisher determined 
that the specific- heat pseudocritical point scales with A = 1/z/ = 1 [32]. However, Ising 
models defined on other two-dimensional lattices with different topologies have shift ex- 
ponents which differ from the inverse correlation length critical exponent (see Ref. [33] for 
a discussion and references therein). This is despite the fact that the critical properties 
on such lattices are the same as for the torus in the thermodynamic limit. Here we ask 
the question: Is there a criterion for A 7^ 1/z/? We find that there is such a criterion - the 
specific- heat amplitude ratio must be 1. 

Although a detailed FSS theory for the impact angle is lacking, for sufficiently large 
L one may expect that is approximated by the angle subtended by the first zero on the 
real axis at the critical point. 

To find the angle <;/>, let L — )■ 00. If, as in most cases, A = l/i/, then tan</) = C, some 
constant value. If A < 1/v then tan</) = 0. But this is impossible as A cannot be less than 
1/v (otherwise v becomes 1/A). 

The only possibility to have A 7^ l/i/, then, is for A > 1/v. In this case tan</) diverges 
as L — )■ 00, so that (f) = i\ j1. From (13. 2p . then, one finds y4_|_ = A_. This is precisely what 
happens in the Ising and Potts models in c? = 2, where the coincidence of the specific-heat 
amplitudes is guaranteed by self-duality. 



Note added in proof: After submission of this work, we became aware of a recent extended 
scaling approach by Campbell and Lundow [M] which yields a universal amplitude ratio 
A^/A_ = 0.540(4). This value is also compatible with the estimate presented herein. 
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